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ABSTRACT 



Context. Turbulent deflagrations of Chandrasekhar mass White Dwarfs are commonly used to model Type la Supernova explosions. In this 

context, rapid rotation of the progenitor star is plausible but has so far been neglected. 

Aims. The aim of this work is to explore the influence of rapid rotation on the deflagration scenario. 

Methods. We use three dimensional hydrodynamical simulations to model turbulent deflagrations ignited within a variety of rapidly rotating 
CO WDs obeying rotation laws suggested by accretion studies. 

Results. We find that rotation has a significant impact on the explosion. The flame develops a strong anisotropy with a preferred direction 
towards the stellar poles, leaving great amounts of unburnt matter along the equatorial plane. 

Conclusions. The large amount of unburnt matter is contrary to observed spectral features of SNe la. Thus, rapid rotation of the progenitor star 
and the deflagration scenario are incompatible in order to explain SNe la. 

Key words. Stars: supernovae: general - Hydrodynamics - Methods: numerical 



1. Introduction 

Driven by the outstanding potential of Type la Supernovae 
(SNe la) to reconstruct the expansion history of the late uni- 
verse, a number of observational campaigns have collected a 
wealth of empirical data about the ligh tcurves, spectra, and host 
populations of the se powerful events dWood-Vasev et al]|2007t 



fective propagation rate is govern ed only by the turbulence pro - 
duced by large-scale instabilities (INiemever & Wooslevll997h . 



Ries s et alj 120071) . This progress is in contrast to the sober- 



ing fact that despite many decades of theoretical and compu- 
tational efforts, we are still lacking a solid understanding of the 
explosion mechanism, or the combination of different mecha- 
nisms, that can explain the class of SNe la in its entirety. While 
the majority of SN la explosions is believed to involve the 
rapid thermonuclear c ombustion of Chandrasekhar mass CO 
White Dwarfs (WDs) faillebrandt & Niemeyerll200ri Uppke 
2008). it is still debated exactly how and where the star ignites 
IWooslev et"aill2004 ISchmidt & Niemeveill20o5 iRopke et aT 



20061) . and whether the combustion front propagates sub 



sonically as a deflagration or supersoni c ally as a detonation 
( Plewa etal] 120041: iRopkeet ail l2007at IRopke & Niemeyei 
2007). However, it is clear that if the front is subsonic, its ef- 
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From a purely theoretical point of view, the turbulent de- 
flagration model is favoured compared with its c ompetitors 
that i nvoke a spontaneous transition to a detonation dNiemeverl 
1 19991) . but large 3D simulations indica te it may be inca - 
pable to explain the full class of SNe la (R6p keetalJl2007ab . 
Nevertheless, because of its conceptual simplicity and lack of 
free parameters, we will use it as a testbed for studying the in- 
fluence of rapid rotation on the outcome of the explosion. Since 
the duration of thermonuclear burning is determined by the dy- 
namical time scale of the WD, the star must rotate nearly crit- 
ically in order for rotation to have any significant impact. This 
appears plausible at first sight, since most progenitor scenarios 
for SNe la involve the accretion of almost a solar mass of ma- 
terial from a binary companion before the thermonuclear run- 
away, allowing the WD to pick up a substantial amount of an- 
gular momentum. Indeed, some evolutionary calculations pre- 
dict critical rotation o f the WD at the time of the explosion 
dYoon & Langerll2005l) . On the other hand, the transport and 
dissipation of angular momentum inside the star are poorly un- 
derstood and direct observations are absent. Hence, we treat 
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the amplitude and shape of the rotation law as essentially un- 
determined and attempt to investigate their influence in a pa- 
rameterised manner. As we will show, our general conclusion 
is robust with respect to the details of the rotation law: nearly 
critical rotation of the progenitor WD together with a pure de- 
flagration model for the combustion front are inconsistent with 
the observations of SNe la. 

The first multidimensional simula tions of rapidly r otatin g 
exploding WDs were carried out by ISteinmetz et al. (1992). 
However, in contrast to our work they concentrated on pure 
detonations instead of deflagrations. We will revisit the prompt 
detonation scenario in a follow-up paper. 

This paper is structured as follows. Section 2 summarises 
the theoretical background on rotation of SNe la progenitor 
stars. The motivation of the employed rotation laws is given. 
Different methods for the numerical initiation of the burning 
process are presented in section 3. Section 4 describes the com- 
putational grid, the flame modelling, treatment of the gravita- 
tional potential, the employed WD models, numerical stability, 
and the parameter study. The results of the numerical study are 
listed in section 5. The impact of rotation on the explosion is 
explained by considering buoyancy effects. Also, the influence 
of model parameters and shear motion caused by rotation is 
presented. An interpretation with respect to spectral features 
finalises the section. Section 6 concludes the paper. 

2. Rotation of CO white dwarfs 

In the c ontext of the single degenerate progenitor scenario for 
SNe la (lLividl2000h . a white dwarf star consisting of carbon 
and oxygen (CO WD) accretes matter from a non-degenerate 
companion via a Keplerian disc. Since the average WD mass 
is ~ 0.6 Mq, the amount of matter that has to be accreted 
prior to C ignition is large (this occurs soonest for non-rotating 
CO WDs close to their Chandrasekhar mass, i.e. at Mwd ~ 
1.4 Mq). In addition to the matter itself, angular momentum 
is accreted. Therefore it is reasonable to assume that the WD 
spi ns up during its accretio n phase. 



Yoon & Langerl (120041) (YL) analyse the physics of accret- 



ing white dwarfs (AWDs) in detail and claim that in case 
of rapid rotation, the limiting mass for C ignition or col- 
lapse can greatly exce ed the value of the Chandrasekhar-mass 
( Yoon & Langerl2005 ). YL's stellar evolution code includes the 
effects of accretion-induced heating and energy transport, an- 
gular momentum transport by various instabilities and the ef- 
fect of rotation on the WD structure by modifications in the 
momentum and energy conservation equations. Electron and 
ion viscosities, relevant for the onset of instabilities, are taken 
into account. Transport of the accreted angular momentum is 
treated as a diffusion process. 

As their main conclusion regarding the angular velocity 
profile of the SN la progenitor, accreting WDs rotate differen- 
tially throughout their evolution. The angular velocity reaches 
a maximum, Q pea k, within the star, increasing from the centre 
outwards but decreasing toward the surface. This behaviour is 
attributed to the fact that the slowly rotating inner part con- 
tracts faster than the rapidly rotating surface layers as the total 
mass of the accreting star increases. Its value is in the range of 



1.5 rad/s < Q pea k 5= 6.5 rad/s when the central density of the 
WD reaches the ignition density of 2.0 x 10 9 g/cm 3 . The cor- 
responding total masses are -1.5 and ~ 2.1 Mq, respectively. 

As YL's results suggest, the rotation law(s) realised by na- 
ture may depend on different factors such as efficiency of ac- 
creted angular momentum transport, the accretion time scale, 
the time scale for loss of angular momentum, the total bi- 
nary mass budget, or magnetic fields. Therefore, in this work 
we study the impact of a parameterised class of rotators with 
different rotation laws and, consequently, different explosion 
masses. Specifically, four types of AWD rotators are consid- 
ered as SN la progenitors. They are - listed in order of in- 
creasing stellar mass - denoted by AWD1, AWD2, AWD4, and 
AWD3 and include the full spread of AWD masses Physical 
quantities and numerical details are itemised in Table Q] The 
central density is p c = 2.0 x 10 9 g/cm 3 for all models. M/„ ? /, 
is the mass for all material if p > 5.248 x 10 7 g/cm 3 , and 
M med if 5.248 x 10 7 g/cm 3 > p > 1.047 x 10 7 g/cm 3 . /? 
denotes the ratio of rotational energy and gravitational bind- 
ing energy, J the angular momentum. In the non-rotating case 
(p c = 2.0 x 10 9 g/cm 3 , M - 1.4 Mq) M high and M med 
are 1.268 Mq (90.6 %) and 1.110 Mq (7.9 %), respectively. 
Furthermore, the non-rotating star possesses the following 
quantities: r = 2.1905 x 10 8 cm, E inl = 2.5228 x 10 51 erg, 
E grm = -3.0275 x 10 51 erg, E mt = 0.000 x 10 50 erg, and 
E b ind = E inl + Eg rao + E rot = -5.047 x 10 50 erg. Since rigid 
rotation is believed to occur in the pre-supernova convective 
core due to very efficient exchange of angular momentum, rigid 
rotation of the interior is assumed for the AWD4 model (see 
Fig. [TJ while the rotators AWD1 to 3 rotate differentially ev- 
erywhere. 

Besides AWD rotation, different rotation laws have been 
tested for pre-supernova WDs. If there are strong braking mag- 
netic torques or efficient viscous angular momentum trans- 
port, a critically rotating rigid rotator might result. In con- 
trast to differentially rotating WDs, the mass of this kind of 
object increases only slightly by 0.07 Mq for a central den- 
sity p c = 2.0 x 10 9 g/cm 3 . In addition, the frequently used 
j - const, and v - const, rotation laws, denoting approximately 
constant specific angular momentum and rotational velocity, 
respectively, are employed (cf. Tabled. The parameter adr de- 
notes the fraction of the equatorial radius up to which rotation 
has approximately constant spec ific angular momentum and 
rotational velocity, respectively (Eriguchi & Mti ller|[l985 ). For 
both of these rotation laws, angular velocity decreases steadily 
from the rotational axis outwards. See section H for details of 
the numerical treatment of the employed rotating WDs. 

Under certain restrictions, it is also possible to interpret 
our model in the context of the double-degenerate scenario 
for SNe la, i.e. the merger and immediate explosion of two 
WDs. Rapid rotation is particularly well motivated in this case. 
However, two conditions need to be satisfied for our model as- 
sumptions to apply: first, in order to reach hydrostatic equi- 
librium of the initial configuration, the explosion must only 
be ignited after a sufficiently long delay of several dynami- 
cal time scales. Second, it is unclear whether in this scenario, 
ignition will occur first near the surface and propagate in- 
ward (see Hillebrandt & Niemeyer (2000) for details and ref- 
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Fig. 1. Density contour lines (left panel) and the rotation law both in spatial (middle panel) and mass coordinates (right panel) for 
the AWD4 model. Unlike the dotted line that denotes the spherical mass coordinate, the solid line line indicates the cylindrical 
mass coordinate in order to describe cylindrical ro tation. Note that the extent of the convective core (i.e., the rigidly rotating inner 
part) is estimated following lYoon & Lan ger (2004). 



erences). Note that while the former is simply a restriction on 
the applicability of our initial conditions, the latter is a serious 
problem for the double-degenerate scenario in itself. For the 
purpose of our study, we remain agnostic about the progeni- 
tor history as long as the conditions for viability of our initial 
model are met. 



3. Ignition conditions 

The explosion is assumed to be a turbulent deflagration 
( Hillebrandt & Niemeyerl2000l) . The structure of the flame sur- 
face at the time of ignition is known to have a strong impact on 
the explosion outcome. In order to investigate the impact of cer- 
tain ignition situations, we used a variety of different methods 
to initiate the burning front in our simulations. 

A standard flame morphology that has be en employed in 



all speed of the process and thus the total number of ignition 
events. 



sever al applications is the c3 configuration dReinecke et al 



20021) . representing a central ignition with three bulges in az- 
imuthal direction per octant. Besides the fact that the c3 ig- 
nition is reproducible and ensures similar conditions even if 
initiated in different initial models, a major advantage of the 
c3 ignition is the applicability for coarse numerical resolut ions. 
As a modification of the static multi-point ignition iRopkd 



(2005), the stochastic ignition model initiates flame seeds 
rando mly at different instants in time dSchmidt & Niemever 
2006). The probability of an ignition depends on the back- 
ground temperature and the magnitude of the convective tem- 
perature fl uctuations which are modelled by a mixing-lengt h 
approach dWoosley et al.1 l2004t IWunsch & Woosleyl Eooi) . 



Thus, the density of the emerging ignitions decreases with the 
distance from the centre and becomes statistically symmetric in 
angular directions. Due to the expansion of the white dwarf, the 
ignition probability decreases with time. The only parameter of 
the model, the exponentiation parameter C e , controls the over- 



An alternative scenario is the dipole jet flow (W oosley et al 



2004), where ignition sites are clustered along the rotation axis 
in one of the two hemispheres. The alignment of ignitions with 
the outward flow direction is explained by the fact that the con- 
dition necessary for ignition is met after stellar material has 
been transported through the hot stellar centre by the dipole jet 
flow. However, the dipole flow - also likely to occur in the non- 
rotating WD as a consequence of the convective pre-supernova 
motion - might be e liminated by a moderate amount of rotation 
dKuhlen et al.ll2006l) . 



4. Numerical Method 



We use the hydrodynamics code Prometheus ( Reineckell2001 ) 
for the simulation of explosions of rotating WDs in three di- 
mensions. The computa tional grid is a moving grid composed 
of two nested sub-grids Ropke (2005). The inner grid captures 
the burned region and the flame surface, while the outer grid 
expands together with the expanding WD. A crucial parameter 
for numerical simulations of SNe la is flame propagation veloc- 
ity, because the evolution of the flame front and, consequently, 
the energy release from thermonu clear burning i s greatl y en- 
hanced by turbulence . Following Schmid t et al. d2006a ) and 
Schmidt et al. ( 2006bl) . the speed of the flame front is computed 
by a subgrid scale model that is based on a balance law for the 
turbulent energy. 

While previous simulations were performed with a gravi- 
tational potential that is spherically symmetric, we employ a 
multipole solver for the computation of the gravitational po- 
tential of rotating stars. A parameter study concerning different 
expansion contributions in rotating WDs showed that <& g {r) is 
approximated with sufficient accuracy if the axisymmetric mul- 
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tipole expansion is truncated after the quadrupole term. The nu- 
meric al implementation w as adapted from the public FLASH 
code dFrvxell et al.ll2000l) . The FLASH implementation of the 
multipole solver is in turn adapted fr om the original implemen- 
tation bv lMiiller & Steinmet3dl995h . 

Rotating density stratifications in hydrostatic equilibrium 
had to be generated as initial conditions for the hydrodynamic 
simulations (see Fig.Q]for an example). We cons t ructed initial 
models using the method of Eriguchi & Mullen (119851) . This 
method, which is an iterative procedure of correcting the den- 
sity field based on local deviations from hydrostatic equilib- 
rium, guarantees the creation of stars in hydrostatic equilib- 
ria for a broad range of rotation laws. The WD equation of 
state was approxim ated using the piecewise polytropic approx- 
imation method (Miille r & Eriguchil ll985). Moreover, the ro- 
tators obey the charac teristic AWD rotation laws as shown by 
Yoon&Langen(l2004 . This was achieved with the implemen- 
tation of a polynomial representation of the angular velocity 
(which was terminated after the 5 th order; cf. Table Q] for the 
chosen coefficients): 



O c for < S < Srigid 

Q c + Ci S + C2 s 2 

+ c 3 s 3 + c A s A + c 5 s 5 for 5 rigid < s < W 

(1) 



Here, s is the normalised distance from the rotation axis and 
Srigid the normalised distance from the rotation axis out to the 
position for which the angular velocity has a constant value of 



equator 



and 



frigid — 



surf 



frigid 
equator 
surf 



(2) 



Since a strong magn etic field might preven t the rapid rotation 
of accreting WD's (lYoon & Lang er 2004), and the possible 
shape of the rotation law may vary to a large degree for differ- 
ent assumptions in the accretion model or simply from star to 
star, the initial models employed for this study were chosen to 
cover a broad range of possible rotation laws and, accordingly, 
masses. Besides the AWD rotation laws, we also considered 
rigid rotation as well as the j - const, and v - const, rotation 
laws, both stating decrea sing angular velocity outwards (cf. 
Eriguchi &Mullerld 1985b for details on both rotation laws). 

Tables Q] and f2] summarise the parameters of our models. 
The masses of the AWD initial models range from 1 .64 Mq to 
2.07 Mq. Model AWD4 contains a rigidly rotating inner core 
(up to the equatorial position s — 6.76 x 10 6 cm) that is thought 
to arise from the pre-supernova convective motion within the 
star. The rotation law for this model is plotted in Fig. Q] Note 
that differential rotators can become heavier than the critical 
rigid rotator since Q. is allowed to have a high value within the 
star and to drop towards the equatorial surface. As a conse- 
quence, overcritical rotation is avoided. 

The rotation law was not fitted for an exact overlap with 
those proposed by YL. In contrast, our parametrisation was 
chosen to emulate the central angular velocity, Q c , the maxi- 
mum of O, and the total mass. Given the uncertainties in the 
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0.401 


M,„, r ,iM r%i 


11.0 
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2.0180 
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1.8833 


1.8836 


r„„„ riO^cml 
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3.2871 


3.2471 


3.3822 


4.1116 


Ei„, riO 51 erel 


2.7395 


2.8895 


2.9148 


2.9604 


E„ r „„ riO 51 erel 


-3.6701 


-4.0330 


-4.1419 


-4.5283 


E nt [10 5U ergl 

*-* rot l o J 


1.748 


2.625 


2.991 


4.465 


E iind [10 iu erg] 


-7.558 


-8.810 


-9.280 


-11.214 




4.7639 


6.5092 


7.2206 


9.8608 


J [10 M, g cm 2 /s] 


0.91512 


1.1968 


1.3602 


2.1951 



Table 1. Physical quantities and numerical coefficients of the 
constructed AWD initial models. 



model 


rigid 


J con st 


V const 


fl c [rad/s] 


0.37242 


6.4076 


10.120 


t equator t ' pol 


1.477 


2.183 


2.183 


adr 




0.5 


0.2 


M[M Q ] 


1.4661 


1.8025 


1.7067 


Mhigh [Mq] 


1.302 


1.534 


1.493 


M high IM [%] 


88.8 


85.1 


87.5 


M med [Mq] 


0.126 


0.191 


0.157 


M med /M [%] 


8.6 


10.6 


9.2 


r po i [10 8 cm] 


2.1721 


1.8087 


1.8134 


r ■equator [10 8 Cm] 


3.2092 


3.9481 


3.95858 


E m [10 51 erg] 


2.5878 


3.0072 


2.9678 


E gmv [10 M erg] 


-3.2074 


-4.2779 


-4.0604 


E m , [10 M) erg] 


0.442 


3.163 


2.393 


E bmi [10 50 erg] 


-5.754 


-9.544 


-8.533 


£[%] 


1.3769 


7.3941 


5.8933 


J [10 5U gcnr/s] 


0.37747 


1.3785 


1.0535 



Table 2. The same parameters as in Table [T] for rigid rotation 
and the "exotic" j const and v consl rotators. 



predictions for rotation laws, this choice is sufficiently realistic 
for our purposes. 

The mass accretion from a Keplerian disc implies a criti- 
cally rotating surface of the AWD rotators which cannot eas- 
ily be incorporated because the initial model algorithm aborts 
when mass shedding for a model is detected. Therefore it is 
assumed that the surfaces of most constructed models in this 
study do not rotate critically. This is not problematic since the 
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- stochastic ignition; here, the ignition speed was varied 
by changing the parameter C e (cf. section[3]) 

All simulations for this work were carried out in three spa- 
tial dimensions. At least one hemisphere or even a full star 
(necessary for the simulation of the dipole jet flow ignition) 
were considered. Since our objective was a parameter study, we 
chose a relatively coarse initial resolution of 7.5 x 10 5 cm on a 
256 2 xl28 grid for one hemisphere. In order to get an estimate 
of the accuracy of all simulations, the "rigid C e — 5 x 10 4 " 
explosion model was redone with doubled resolution. It was in 
remarkable agreement with the run on a normal grid (the ex- 
plosion energetics differ by less than 4%, the produced species 
differ by less than 2%). This result indicates that the resolu- 
tion employed for most of the simulations is sufficient since 
the explosion energetics and species composition are almost 
identical. 



Fig. 3. Gravitational acceleration g gmv (dotted line) and the ef- 
fective gravitational acceleration gr e ff (solid line) along the ro- 
tation axis, the stellar diagonal, and the equatorial plane for the 
AWD4 rotator. The edges of the red line (at ~ 3000 km) and the 
blue line (at ~ 3400 km) appear due to the absence of centrifu- 
gal acceleration outside the star. The dashed thick black line 
represents the spherically symmetric gravitational acceleration 
for the non-rotating star of equal central density. 



5. Results 

In the following, we will first analyse deflagration for the 
AWD4 rotator, which is ignited by stochastic ignition. The 
strong impact of rotation on the flame morphology can be re- 
garded as typical for the deflagration of rapidly rotating WDs. 
Variations of the stochastic ignition process and the influence 
of the rotation laws are investigated next. We will also comment 
on the role of shear motion and the expected spectral features. 



outer layers of a star do not affect the nuclear burning in the 
interior. 

We tested whether the initial models remain hydrodynam- 
ically stable on the grid structure used for the SNe la simula- 
tions. Hydrodynamical simulations were followed to t — 1.5 s, 
the time when burning would have stopped if the WD had 
been ignited. A comparison of density contours at different in- 
stants served as a simple check for hydrostatic stability. Up to 
t — 1 .0 s, the contour lines retained their shape except at the 
surface. At later times, the WD's core density also decreases 
somewhat because of the limited feasibility to perform circu- 
lar motion on a rectangular cartesian grid. Two reasons are 
responsible for amplified deviations in the outer layers. First, 
the resolution of the grid is lower outside the core region and 
second, the pressure gradient is neglected below the threshold 
density of p — 10~ 3 g/cm 3 . A more quantitative test was done 
by looking at the radial velocities emerging during the test of 
hydrostatic stability. Radial velocities in the inner region are at 
all instants found to be < 4 x 10 7 cm/s and therefore one order 
of magnitude lower than the characteristic velocities after ig- 
nition. In conclusion, stability of all discretised models on the 
grid can be assumed. 

In order to cover a broad range of possible explosion sce- 
narios, we varied the following parameters: 

- the initial model: norot, rigid, AWD1, AWD2, AWD4, 
AWD3, jconst, and u cons , 

- the ignition scenario: 

- static ignition (c3) 



5.1. Anisotropic burning induced by rotation 

The existence of a preferred direction for the flame propaga- 
tion in rotating white dwarfs can be attributed to two effects. 
First, the ashes of SNe la are subject to buoyant motion, which 
is stronger in polar direction. Second, for convective motion 
in equatorial direction, the burning products have to gain an- 
gular momentum, whereas the opposite is true for inflowing 
fuel. This effectively gives rise to an angular momentum bar- 
rier, which suppresses the turbulent spreading of the deflagra- 
tion front in the direction orthogonal to the rotation axis. 

Fig. |2(a)| shows the gravitational acceleration g gmv inside 
the AWD4 rotator. As a result of the larger density gradient 
along the rotation axis, gr arav is stronger (indicated by longer 
arrays) in this direction compared to the equatorial plane. 
Buoyancy is caused by the effective gravitational acceleration 
g e s, i.e., the gravitational acceleration including the centrifu- 
gal acceleration as plotted in Fig. |2(b)| g e s is orthogonal to the 
coinciding isobaric and isopycnic surfaces and exactly com- 
pensates the pressure gradient in hydrostatic equilibrium. The 
difference in g e s along the polar and equatorial direction is even 
larger than the difference in g mv , which can be seen in Fig. [3] 

For the veloci ty Vb of a buoyant bubble, the following 
expression holds ( Davies & Tavloi 1950b Glimm & Li 1988; 
Readlll984h : 



(3) 



where the dimensionless constant C\ is about 0.5, D is the bub- 
ble diameter, and the Atwood number At = (pf -pb)/(pf + Pb) 
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Fig. 2. Density contour lines together with the gravitational acceleration gr grav (a) and the effective gravitational acceleration g e g 
(b) that is indicated by arrows, respectively, for the AWD4 rotator. Note that gr e ff is orthogonal to the coinciding isobaric and 
isopycnic surfaces (outside the star, the centrifugal acceleration is zero). 



specifies the density contrast between the rising burning prod- 
ucts of density pb in fuel of density pf > pt,. The dependence of 
Vb on the position inside the star is mainly given by the varia- 
tion of g e ff . The magnitude of g e g differs by a factor up to 1 .3 in 
polar and equatorial directions, (see Fig. |3). From equation (O 
follows that the rising velocity of ashes in polar direction can 
exceed the velocity in equatorial direction by 15 %. The ve- 
locity of the fastest bubbles is of the order ~ 10 8 cm/s (with 
At = 0.5, Igfeffl = 2 x 10 10 cm/s 2 , and a bubble diameter of 
10 s cm), a value that is approximately 10 % of the speed of 
sound. Accordingly, smaller bubbles rise more slowly. A lower 
limit for the rising speed is of the order of 10 6 cm/s (with 
At = 0.1, a lower gravitational acceleration and the bubble di- 
mension of 10 5 cm). 

In the following, we consider the explosion of the AWD4 
rotator in the stochastic ignition scenario with C e — 5 x 10 4 . 
Fig.|4]shows the flame surface (represented by thick white con- 
tours) embedded in density profiles at four instants. The simu- 
lation was carried out for one hemisphere only. Thus the con- 
tour plots — cross sections in the i/z-plane along the rotation 
axis — exhibit equatorial symmetry. During the first 0.1 sec- 
onds, 363 ignition spots are set isotropically in the interior of 
the star. Even at this early period, an agglomeration of the flame 
around the rotation axis (which points along the x axis) is vis- 
ible. Subsequently, the anisotropy of the flame grows because 
of the effects discussed above. In addition, the turbulent flame 
propagation speed is amplified in the polar direction due to the 
stronger production of turbulence. The thin contour lines (red 
in the online version) in Fig. |4] enclose densities greater than 
5.248 x 10 7 g/cm 3 (threshold for producing iron group ele- 
ments) and 1.047 x 10 7 g/cm 3 (threshold for producing inter- 
mediate mass elements). 



In conclusion, both the buoyancy that is enhanced at the 
poles in rotating WDs and the blocked mixing orthogonal to 
the rotation axis lead to a strongly anisotropic flame that leaves 
behind unburnt material at the stellar centre. 



5.2. Ignition speed in the stochastic ignition scenario 

As outlined in section [3] the stochastic ignition process con- 
tains one free parameter: the exponentiation parameter C e that 
controls the overall speed of the creation of ignit ions. Just as in 
the original work by Schmidt & Niemever ( 20061) . this parame- 
ter was varied over a broad range in order to study the influence 
of C e . For this reason, the AWD3 rotator was ignited for C e be- 
tween 5x10' and 5 x 10 7 (see Table [3j> in the isotropic mode, 
i.e., for ignitions that are placed without preference. Q denotes 
the total spherical angle covered by the simulation domain. 

As summarised in Table [3] increasing C e results in a 
higher amount of iron group elements (IGEs), while the 
production of intermediate mass elements (IMEs) does not 
ch ange significantly. However, similar to the results obtained 
by lSchmidt & Niemeyerl l 20061) . the energy output and the IGE 
production do not increase further for C e > 5 x 10 . This effect 
can be explained by the saturation or even decrease of the flame 
surface during the ignition phase for an extrem ely higher num- 
ber of ignition events (see lRopke et al.l (120061) for a discussion 
of the same trend for the static multi-spot ignition). However, 
employing the value C e — 5 x 10 5 amounts to burning the en- 
tire core almost instantaneously because within the first tenth 
of a second almost 10 3 ignitions are initiated within a central 
region of the typical radius of ~ 350 km. Since high values of 
C e not only seem unphysical but also cannot serve as a remedy 
for the problems arising for the stochastic ignition in combina- 
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Fig. 4. Density contour lines for the pure deflagration of the AWD4 rotator at different instants for stochastic ignition (C e = 
5 x 10 4 ). Cross sections along the rotation axis of these simulations exhibiting equatorial symmetry are shown. 



c e 


5 x 10 1 


5 x 10 3 


5x 10 4 


5 x 10 5 


5 x 10 7 


10 4 


mode 


iso 


iso 


iso 


iso 


iso 


dipole 


Q [n] 


2 


2 


2 


2 


2 


4 


hn 


13 


152 


391 


809 


19980 


247 


t = 5s 


E km [10 iu erg] 


2.87 


3.61 


5.30 


6.90 


6.11 


4.43 


E m [10 50 erg] 


0.21 


1.33 


3.54 


5.47 


4.62 


2.55 


E nuc [10 51 erg] 


1.050 


1.182 


1.420 


1.619 


1.535 


1.310 


IGEs [MO] 


0.40 


0.51 


0.69 


0.79 


0.74 


0.62 




20 


25 


34 


39 


37 


31 


IMEs [MO] 


0.51 


0.46 


0.41 


0.46 


0.45 


0.40 


[%M tot ] 


25 


23 


20 


23 


23 


20 


C+O [MO] 


1.10 


1.04 


0.91 


0.76 


0.82 


0.99 


[%M m ] 


55 


52 


46 


38 


40 


49 



Table 3. Energetics and compositions for the deflagration mod- 
els activated by stochastic ignition with different values for C e 
in the AWD3 rotator (M tol = 2.02 Mq). 



tion with rotation (cf. section |5~.3.2| i, the values of C e — 5 x 10 3 
and C e — 5 x 10 4 were used in this work. 



The dipole jet flow scenario was also tested by means of 
the AWD3 rotator. However, as the impact of buoyancy al- 
ready suggests (cf. section 15.11 ). this ignition realisation does 
not ameliorate the problem of remaining fuel close to the cen- 
tre. If ignitions are set mainly in one hemisphere (without loss 
of generality, the southern one), the northern hemisphere is also 
burnt at a later time, but C and O are left over in the core. In 
the non-rotating case, a strong off-centre ignition also leads to 
a delayed burning of the distant stellar p art, but less centrally 
located fuel remains (Ropkeetal ] l2007bl) . 



We also found that the detailed location of randomly cho- 
sen ignition points has a significant impact on the outcome of 
the simulation. The "rigid C e — 5 x 10 3 " explosion model was 
simulated twice and is summarised in the two colums rigid and 
rigid * of Table [5] The only difference was a different seed for 
the random number generator which sets the ignition position 
within a shell whose number is determined by the Poisson pro- 
cess. Due to the randomness of the ignition process, the total 
number of ignitions is not entirely but approximately the same: 
167 vs. 163 ignitions were set during the simulations. However, 
since bubbles inserted close to the rotation axis grow much 
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norot 


rigid 


AWD1 


AWD2 


AWD4 


AWD3 


M M [MO] 


1.40 


1.47 


1.64 


1.74 


1.79 


2.02 


Q. [n] 


2 


2 


2 


2 


2 


2 


t = 0s 


Egrav [10 M erg] 


-3.027 


-3.204 


-3.665 


-4.027 


-4.135 


-4.518 


E ro , [l(P u erg] 


0.00 


0.44 


1.75 


2.62 


2.99 


4.46 


yS[%] 


0.00 


1.38 


4.78 


6.50 


7.22 


9.87 


t = 5s 


E km [10 5U erg] 


6.49 


6.53 


6.63 


5.99 


5.42 


5.28 


E m [10 M, erg] 


5.89 


5.87 


5.74 


4.84 


4.09 


3.50 


E mlc [10 31 erg] 


1.076 


1.129 


1.291 


1.318 


1.289 


1.419 


IGEs [MO] 
[%M tot ] 


0.54 

38 


0.56 
39 


0.65 
40 


0.65 
37 


0.64 
35 


0.67 

33 


IMEs [MO] 
[%M m ] 


0.28 
20 


0.30 
20 


0.34 
20 


0.36 
21 


0.35 
20 


0.44 

22 


C+O [MO] 
[%M tot ] 


0.58 
42 


0.61 
41 


0.65 
40 


0.73 
42 


0.80 
45 


0.90 
45 



Table 4. Energetics and compositions for the deflagration 
model with c3 ignition in different rotators. 



more rapidly and tend to dominate the explosion, the gener- 
ated nuclear energy differs by ~ 20 % (the more fuel remains 
unburnt, the lower is the value of produced nuclear energy). 
This demonstrates the high sensitivity of the rotating WD de- 
flagration scenario on initial conditions. The sensitivity to the 
initial distribution of the seed flame was also observed for non- 
rotating WDs but in that con text it caused only small deviations 
( Schmidt & Niemeyerl2006 ). Any stellar rotation enhances this 
finding as a result of the oblate shape of the rotators. 



5.3. Variation of the rotation law 

Differences in the rotation law give rise to substantial variations 
of the total mass of the thermonuclear supernova progenitor. In 
particular, rapidly rotating WDs can greatly exceed the canon- 
ical Chandrasekhar mass: the progenitor stars under considera- 
tion possess masses in the range of 1 .40 Mq (no rotation) and 
2.02 Mq (model AWD3). An important conclusion can already 
be drawn by looking at the mass densities which the different 
WD models exhibit. The total mass increases by ~ 44 % from 
the non-rotating star to the AWD3 rotator. However, the mass 
fraction with densities greater than 5.25 x 10 7 g/cm 3 , for which 
IGEs are produced, grows by only ~ 16 %. The reason is that 
the fraction of dense material within the star decreases from 
~ 91 % in the non-rotating star to only ~ 73 % in the AWD3 ro- 
tator which — although heavier by ~ 0.6 Mq — contains only 
~ 0.2 Mq more IGEs. This is illustrated in Figs. [5(a)1 and [5(b)l 
Therefore it is doubtful whether the pure deflagration of rapidly 
rotating super-Chandrasekhar-mass models yields significantly 
more 56 Ni compared to the non-rotating case. We consider the 
different rotation laws for the c3 and stochastic ignition scenar- 
ios in turn. 




-10 - 



t[s] 



(a) rotation law: c3 ignition 




(b) rotation law: stochastic ignition 

Fig. 7. Temporal evolution of the total energy for the deflagration 
study that is initiated by c3 ignition in different rotators (a), and 
by stochastic ignition in different rotators (C e — 5 x 10 4 for the 
AWD4 rotator, else C e = 5 x 10 3 ) (b) 



5.3.1 . c3 ignition 



Fig. |7(a)| shows the temporal evolution of the total energy E t0 , 
for different rotators that are ignited with the c3 ignition. Note 
that the more massive WDs start from more negative values of 
total energy. This is caused by the larger gravitational energies 
Egrav for heavier stars (see Table HJ. All simulations end in ex- 
plosions that cause the star to unbind at t ~ Is. This takes 
longest for the heaviest WD, AWD3. The differences in E tot 
that are apparent in the beginning become smaller during the 
explosion process, which is caused by a higher amount of re- 
leased nuclear energy E mlc for the more massive WDs. 

The non-rotating star exhibits the highest value of E m in 
the homologous expansion phase. E m at t — 5 s consists mainly 
of kinetic energy Ekm of the ejecta for the typical temporal evo- 
lution of the energy contributions. Accordingly, the ejecta are 
fastest for the non-rotating progenitor on average. 

With respect to the stellar composition after the explo- 
sion, the following trend is visible from Tableland Fig. |6(b)| 
the heavier the progenitor, the more fuel is left unburnt. The 
amount of IGEs after 5 s ranges from 0.54 Mq for the non- 
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1,4 1,5 1,6 1,7 1,8 1,9 2 1,4 1,5 1.6 1,7 1,8 1,9 2 



M w D [ M «J M wd[ M .J 

(a) no expansion: absolute fractions (b) no expansion: relative fractions 

Fig. 5. Estimate of the burning products if the stars were not expanding during the burning process. The figure shows the fractions 
of the total mass in the high density regime (p > 5.25 x 10 7 g/cm 3 , "IGEs"), the medium density regime (5.25 x 10 7 g/cm 3 >p> 
1.05 X 10 7 g/cm 3 , "IMEs"), and the low density regime (p < 1.05 x 10 7 g/cm 3 , "C+O"). 




Fig. 6. Burning products after t — 5 s as the outcome of simulations that are ignited by stochastic ignition (solid lines) and 
c3 ignition (dashed lines). Note the influence of the stochastic burning on the explosion outcome as can be seen for the rigid 
rotator: there the burning products are taken as the mean of two identical simulations that differ only in the random numbers for 
stochastic ignition. 



rotating star to 0.67 Mq for the rapidly rotating model AWD3. 
However, the fraction of IGEs with respect to the total mass 
decreases for rapid rotation. Only 33 % of the total 2.02 Mq is 
converted to IGEs for the AWD3 rotator. The amount of IMEs 
rises for more rapid rotation, its fraction is remarkably constant 
at ~ 20 % for all rotation laws. Combined with the fact that the 
amount of unburnt material becomes larger for faster rotators, 
this result indicates that the higher absolute amount of IMEs for 
the rapid rotators is not obtained "at the cost" of fuel in contrast 
to observations. It is a consequence of the small amount of pro- 
duced IGEs. 



5.3.2. Stochastic ignition 

Basically, the trends found in the case of the c3 ignition 
scenario remain valid for explosions with stochastic ignition. 
However, the energetics of the explosions show a larger vari- 
ation (see Fig. |7(b)| and Table [5]), and even less IGEs are pro- 
duced while more unburnt fuel is left over (see Figs. |6(a)| and 
|6(b)| as well as Table [5]) . This is explained by the fact that, un- 
like the symmetric c3 ignition, the stochastically ignited bub- 
bles are immediately subject to buoyant motion since the net 
buoyancy is not balanced by a symmetric alignment of the 
ashes, therefore having less time to burn the central region ef- 
ficiently. Also fluctuations in the composition between differ- 
ent progenitor models are more distinct for stochastic ignition 
compared to c3 ignition. This is related to the fact that not only 
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the overall number of ignitions but also the specific locations 
of the ignitions influence the explosion outcome. 

The most energetic explosions with high fractions of IGEs 
result from the non-rotating and the AWD1 progenitor models. 
For the latter model, the central angular velocity is low com- 
pared to the other differentially rotating models. Since convec- 
tion in the pre-supernova phases causes efficient angular mo- 
mentum transport throughout the core, however, the AWD1 
rotator appears to be a rather unlikely outcome of an evolu- 
tionary scenario leading to a thermonuclear supernova. For the 
"AWD4 C e = 5 x 10 4 " explosion model, on the other hand, 
the amount of IMEs is somewhat higher, but this is compen- 
sated by little production of IGEs. Consequently, a larger frac- 
tion of fuel remains unburnt. The "j const C e = 5 x 10 3 " and 
"Vconst C e = 5 x 10 3 " explosion models generate yet smaller 
amounts of IGEs (see Tabled. However, neither the j co „st nor 
the Vconst rotation law are likely to occur in WDs. For the crit- 
ically rotating rigid rotator, on the other hand, approximately 
the same results are obtained as in the case of the non-rotating 
star. Consequently, this type of rotating WD appears to be a 
progenitor candidate for type la supernovae caused by pure de- 
flagrations. 





J cons! 


» -const 


M m [MO] 


1.80 


1.71 


C e 


5 x 10 3 


5 x 10 3 


n [n] 


2 


2 


h, 


159 


158 


t = 0s 


E grm [10 M erg] 


-4.271 


-4.056 


Erot [10 3U erg] 


3.16 


2.39 


pm 


7.40 


5.89 


t = 5 s 


E km [10 M, erg] 


1.66 


1.29 


Em [10 5U erg] 


-0.75 


-1.11 


E, mc [10 51 erg] 


0.817 


0.679 


IGEs [MO] 


0.34 


0.29 


[%M M ] 


19 


17 


IMEs [MO] 


0.34 


0.27 


[%M M ] 


19 


16 


C+O [MO] 


1.12 


1.15 


[%M tot ] 


62 


67 



Table 6. Energetics and compositions for the deflagration 
model activated by stochastic ignition with C e — 5 x 10 3 in 
the j const and v CO nst rotators. 



5.4. Shear motion during the explosion 

It is a common belief that the shear introduced by differential 
rotation can enhance the explosion strength by increasing the 
flame surface. Although rotating WDs exhibit high rotation ve- 
locities (e.g., the fastest regions of the AWD1 rotator move with 
10 % of the speed of light), within the entire span of time during 
which burning takes place (t < 1.5 s) the direct influence of ro- 
tation on the flame surface is limited. The rotators accomplish 
only about half a rotation if the burning induced expansion is 
neglected. The latter arranges a slowdown of the star as soon 
as the thermonuclear runaway is initiated as a result of angular 
momentum conservation. 

Unless a distinct jump in Q occurs — which is not the case 
for smooth differential rotation — , the flame will not be signif- 
icantly influenced by rotation in a direct way (but, as described 
in section IBTT1 by buoyancy and angular momentum effects). 



Yoon & Langerl (120041) proposed that a rotation could trig- 
ger a deflagration-to-detonation transition since the outer lay- 
ers still rotate rapidly when the outward burning flame causes a 
slowdown of the WD. The resulting gradient in velocity due to 
burning is higher than by differential rotation and could pos- 
sibly generate a detonation if the flame surface is increased 
abruptly. In particular, such a rapid jump of angular velocity i s 
found for the rotation law of a WD merger dYoon et alj|2007l) . 
We do not consider DDTs in this work; however, an investi- 
gation of shear-induced DDT s with the method introduced by 
iGolombek & Niemeverl J2005b would be possible and interest- 
ing. In particular, if observational evidence for rapid rotation of 
the progenitor star should be found, the DDT mechanism will 
be favored since prom pt detonations of ra pid rotators cannot 
explain normal SNe la (Pfann es et al ] l2009h . 



5.5. Expected spectral features 

An inspection of the composition and the kinematics of the 
ejecta in the homologous expansion phase allows a prediction 
of spectral features even without employing a detailed postpro- 
cessing study. Regarding rotation, the most serious discrepancy 
between the simulation results and observations is the predic- 
tion of C and O at low radial expansion velocities, i.e., close to 
the centre. This result, already a problem for the non-rotating 
deflagration model, grows more acute for rapid rotation. Fig. [8] 
shows the total and fractional mass densities after t — 5 s for 
the "norot c3" and 'AWD3 c3" explosion models. For the sake 
of demonstration, only the c3 ignition is presented here. The re- 
sults are qualitatively in accordance with those of the stochastic 
ignition. 

SN la observations indicate that the IMEs should be located 
in the outer regions and enclose the IGEs for normal SNe la, 
whereas C and O must not be present near the centre at all. 
However, as a result of the turbulent combustion whose driving 
force is the Rayleigh-Taylor instability, the species are mixed 
throughout the stellar interior. Furthermore, the existence of a 
significant amount of unburnt fuel in the core is inconsistent 
with observations. The following trend that becomes stronger 
for more rapid rotation emerges: Initially, burning mostly ad- 
vances towards the stellar poles and generates IGEs. Later, 
when the flame spreads toward the equatorial plane and the 
density has already dropped, IMEs are produced. However, 
the burning process ceases before the central part and mate- 
rial close to the equatorial plane can be burnt. This is reflected 
by the distinct peak of C and O at low radial velocities that 
can be seen in Fig. [9] for the AWD4 rotator. For comparison, 
the distributions in radial velocity space are also plotted for the 
non-rotating stochastic ignition model. 
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norot 


rigid / rigid* 


AWD1 


AWD2 


AWD4 


AWD3 


M m [MO] 


1.40 


1.47 


1.64 


1.74 


1.79 


2.02 


c e 


5 x 10 3 


5 x 10 3 


5 x 10 3 


5 x 10 3 


5x 10 4 


5 x 10 3 


Q [n] 


2 


2 


2 


2 


2 


2 


I* 


172 


167/163 


168 


154 


404 


152 


t = 0s 


E gnw [10 M erg] 


-3.027 


-3.204 


-3.665 


-4.027 


-4.135 


-4.518 


E m , [10 5U erg] 


0.00 


0.44 


1.75 


2.62 


2.99 


4.46 


0[%] 


0.00 


1.38 


4.78 


6.50 


7.22 


9.87 


t = 5 s 


E m [10 il, erg] 


8.69 


7.70/5.63 


8.64 


3.67 


3.68 


3.61 


E m [10 M) erg] 


8.16 


7.04/4.82 


7.84 


1.95 


1.97 


1.33 


E mlc [10 51 erg] 


1.282 


1.234/ 1.005 


1.495 


1.016 


1.073 


1.182 


IGEs [MO] 
[%M m ] 


0.67 
48 


0.63/0.50 
43/34 


0.79 
48 


0.47 
27 


0.52 
29 


0.51 
25 


IMEs [MO] 
[%M tot \ 


0.28 
20 


0.30/0.27 
20/19 


0.32 
20 


0.34 
20 


0.31 
17 


0.46 
23 


C+O [MO] 
[%M, o; ] 


0.45 

32 


0.54 / 0.70 
37/47 


0.53 
32 


0.93 
53 


0.96 
54 


1.04 

52 



Table 5. Energetics and compositions for the deflagration model activated by the stochastic ignition with C e 
5 x 10 4 for AWD4) in different rotators. 



5 x 10 3 (C e 



Mass density (t = 5.0 s) 



Mass density (t = 5.0 s) 



i • 



total 

7T 



JVC 



o z 



oo 



20000 4O000 



(a) "norot c3"explosion model 



(b) 'AWD3 c3" explosion model 



Fig. 8. Total (upper left) and fractional (upper right: IGEs; lower left: unburnt fuel; lower right: IMEs) mass densities after 
t = 5 s for the "norot c3" (a) and "AWD3 c3" (b) explosion models. A mixing of the species throughout the star emerges in 
both scenarios, particularly the rapid rotator leaves fuel in the stellar core. Cross sections along the rotation axis of simulations 
exhibiting equatorial symmetry are shown. 



The highest radial velocities in the non-rotating case are 
< 15 x 10 8 cm/s. Consequently, high velocity features can- 
not be explained by the deflagration model, although higher 
velocities might become apparent with increasing resolution. 
In the rotating case, the radial velocities are even smaller (< 
10 x 10 8 cm/s), despite the fact that rotating progenitors in- 
clude a large amount of rotational kinetic energy. This is due 
to the greater gravitational attraction of the rotating stars (cf. 



Table |3), and in addition the released energy is comparatively 
small because of incomplete burning. Fig.[l0]compares the ini- 
tially available rotational energy (bottom of the bars) and the 
kinetic energies of the ejecta after t = 5 s (top of the bars) 
for different WD models. Except for the AWD1 rotator without 
convective core, the disposable kinetic energies become suc- 
cessively smaller for increasing strength of rotation. Note that 
for the "AWD3 C e - 5 x 10 3 " explosion model, the kinetic 
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energy after t — 5 s is even less than the kinetic energy initially 
present in the rotation motion. 



3.0-10" 
2.5-10 2J 



iron group elements (IGEs) 
intermediate mass elements (IMEs) 
carbon and oxygen (C+O) 




2.0-10 8 4.0-10" 6.0-10" 8.0-10" 1.0-10" 1.2-10" 1.4-10" 

v,. [cm/s] 

(a) "norot C e = 5 x 10 4 " 



iron group elements (IGEs) 
intermediate mass elements (IMEs) 
carbon and oxygen (C+O) 




2.0-10° 4.0-10* 6.0-10" 8.0-10" 1.0-10" 1.2-10" 1.4-10" 

v,. [cm/s] 

(b) "AWD4 C e = 5 x 10 4 " 

Fig. 9. Probability density functions in radial velocity space for 
the "norot C e = 5 x 10 4 " (a) and "AWD4 C e = 5 x 10 4 " (b) 
explosion models 



6. Conclusions 

In this paper, we investigated the influence of nearly critical ro- 
tation of the progenitor star on SN la explosions. We performed 
three dimensional numerical simulations of thermonuclear de- 
flagrations for different rotation laws of the white dwarf pro- 
genitors and various initial conditions. The influence of the ig- 
nition process turns out to be larger in rotating stars compared 
to their non-rotating counterparts, which might contribute to 
the observed diversity among SNe la. 

The main result of this study is that the amount of iron 
group elements is not significantly increased compared to the 
non-rotating case, although the rotating WD models are no- 
tably heavier than their non-rotating counterparts. Due to the 
centrifugal expansion, rotators contain more material at low 




Srigic] 



1UAWD1 11AWD2 14AWD4 20AWD3 



Fig. 10. Initially available rotational energy (bottom of the 
bars) and kinetic energy of the ejecta after t = 5 s (top of the 
bars) for successively heavier WD models 



densities in general, thus the amount of material capable of 
burning to iron group elements is only weakly increased if the 
star is burnt by a deflagration. Furthermore, due to rotationally 
induced anisotropic buoyancy effects and, at the same time, 
inhibition of effective mixing parallel to the equatorial plane, 
the flame preferentially propagates toward the stellar poles. 
This leads to comparably weak and anisotropic explosions that 
leave behind unburnt material at the centre and in the equatorial 
plane. Moreover, no significant effects of shear motion acting 
on the flame are observed. Deflagrations of rapid rotators do not 
cause high velocity features but result in overall low expansion 
velocities. As a consequence, the pure deflagration scenario is 
ruled out in the case of rapid rotation of the progenitor star. The 
incineration of the critical rigid rotator shows similar features 
as the non-rotating scenario. Therefore, deflagrations could be 
possible for critical rigid rotation. 

In conclusion, rotation of the progenitor star is unlikely to 
be the parameter that causes the observed variation in peak lu- 
minosities among SNe la. Rapid rotation of the progenitor star 
derived by an accretion study can lead to a great variety in ex- 
plosion strengths in the deflagration scenario, but only within 
a range that is ruled out for observational reasons. Otherwise, 
even the possible critical rigid rotation can not account for a 
significant spread in the explosion outcome. 
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